Search for: All records

Globular clusters (GCs) are particularly efficient at forming millisecond pulsars. Among these pulsars, about half lack a companion star, a significantly higher fraction than in the Galactic field. This fraction increases further in some of the densest GCs, especially those that have undergone core collapse, suggesting that dynamical interaction processes play a key role. For the first time, we createN-body models that reproduce the ratio of single-to-binary pulsars in Milky Way–like GCs. We focus especially on NGC 6752, a typical core-collapsed cluster with many observed millisecond pulsars. Previous studies suggested that an increased rate of neutron star binary disruption in the densest clusters could explain the overabundance of single pulsars in these systems. Here, we demonstrate that binary disruption is ineffective and instead we propose that two additional dynamical processes play dominant roles: (1) tidal disruption of main-sequence stars by neutron stars and (2) gravitational collapse of heavy white dwarf binary merger remnants. Neutron stars formed through these processes may also be associated with fast radio bursts similar to those observed recently in an extragalactic GC.

 

We performed radio searches for the “spider” millisecond pulsar (MSP) candidates 4FGL J0935.3+0901, 4FGL J1627.7+3219, and 4FGL J2212.4+0708 using the Green Bank Telescope in an attempt to detect the proposed radio counterpart of the multi-wavelength variability seen in each system. We observed using the VEGAS spectrometer, centered predominantly at 2165 MHz; however, we were also granted observations at 820 MHz for 4FGL J1627.7+3219. We performed acceleration searches on each data set usingPRESTOas well as additional jerk searches of select observations. We see no evidence of a radio counterpart in any of the observations for each of the three systems at this time. Additional observations, perhaps at different orbital phases (e.g., inferior conjunction), may yield detections of an MSP in the future. Therefore, we urge continued monitoring of these systems to fully characterize the radio nature, however faint or variable, of each system.

 

The Fourier domain acceleration search (FDAS) is an effective technique for detecting faint binary pulsars in large radio astronomy data sets. This paper quantifies the sensitivity impact of reducing numerical precision in the graphics processing unit (GPU)-accelerated FDAS pipeline of the AstroAccelerate (AA) software package. The prior implementation used IEEE-754 single-precision in the entire binary pulsar detection pipeline, spending a large fraction of the runtime computing GPU-accelerated fast Fourier transforms. AA has been modified to use bfloat16 (and IEEE-754 double-precision to provide a “gold standard” comparison) within the Fourier domain convolution section of the FDAS routine. Approximately 20,000 synthetic pulsar filterbank files representing binary pulsars were generated using SIGPROC with a range of physical parameters. They have been processed using bfloat16, single-precision, and double-precision convolutions. All bfloat16 peaks are within 3% of the predicted signal-to-noise ratio of their corresponding single-precision peaks. Of 14,971 “bright” single-precision fundamental peaks above a power of 44.982 (our experimentally measured highest noise value), 14,602 (97.53%) have a peak in the same acceleration and frequency bin in the bfloat16 output plane, while in the remaining 369 the nearest peak is located in the adjacent acceleration bin. There is no bin drift measured between the single- and double-precision results. The bfloat16 version of FDAS achieves a speedup of approximately 1.6× compared to single-precision. A comparison between AA and the PRESTO software package is presented using observations collected with the GMRT of PSR J1544+4937, a 2.16 ms black widow pulsar in a 2.8 hr compact orbit.

 

Abstract Close encounters between neutron stars and main-sequence stars occur in globular clusters and may lead to various outcomes. Here we study encounters resulting in the tidal disruption of the star. Using N -body models, we predict the typical stellar masses in these disruptions and the dependence of the event rate on the host cluster properties. We find that tidal disruption events occur most frequently in core-collapsed globular clusters and that roughly 25% of the disrupted stars are merger products (i.e., blue straggler stars). Using hydrodynamic simulations, we model the tidal disruptions themselves (over timescales of days) to determine the mass bound to the neutron star and the properties of the accretion disks formed. In general, we find roughly 80%–90% of the initial stellar mass becomes bound to the neutron star following disruption. Additionally, we find that neutron stars receive impulsive kicks of up to about 20 km s −1 as a result of the asymmetry of unbound ejecta; these kicks place these neutron stars on elongated orbits within their host cluster, with apocenter distances well outside the cluster core. Finally, we model the evolution of the (hypercritical) accretion disks on longer timescales (days to years after disruption) to estimate the accretion rate onto the neutron stars and accompanying spin-up. As long as ≳1% of the bound mass accretes onto the neutron star, millisecond spin periods can be attained. We argue the growing numbers of isolated millisecond pulsars observed in globular clusters may have formed, at least in part, through this mechanism. In the case of significant mass growth, some of these neutron stars may collapse to form low-mass (≲3 M ⊙ ) black holes. 

The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between 2021 April 16 and 17 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multifrequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment and the 100 m Green Bank Telescope in a 3 yr period encompassing the shape change event, between 2020 February and 2023 February. As of 2023 February, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying time-of-arrival residuals display a strong nonmonotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency,ν) nor a change in dispersion measure alone (which would produce a delay proportional toν⁻²). However, it does bear some resemblance to the two previous “chromatic timing events” observed in J1713+0747, as well as to a similar event observed in PSR J1643−1224 in 2015.

 

Context.By providing information about the location of scattering material along the line of sight (LoS) to pulsars, scintillation arcs are a powerful tool for exploring the distribution of ionized material in the interstellar medium (ISM). Here, we present observations that probe the ionized ISM on scales of ∼0.001–30 au.Aims.We have surveyed pulsars for scintillation arcs in a relatively unbiased sample with DM < 100 pc cm⁻³. We present multifrequency observations of 22 low to moderate DM pulsars. Many of the 54 observations were also observed at another frequency within a few days.Methods.For all observations, we present dynamic spectra, autocorrelation functions, and secondary spectra. We analyze these data products to obtain scintillation bandwidths, pulse broadening times, and arc curvatures.Results.We detect definite or probable scintillation arcs in 19 of the 22 pulsars and 34 of the 54 observations, showing that scintillation arcs are a prevalent phenomenon. The arcs are better defined in low DM pulsars. We show that well-defined arcs do not directly imply anisotropy of scattering. Only the presence of reverse arclets and a deep valley along the delay axis, which occurs in about 20% of the pulsars in the sample, indicates substantial anisotropy of scattering.Conclusions.The survey demonstrates substantial patchiness of the ionized ISM on both astronomical-unit-size scales transverse to the LoS and on ∼100 pc scales along it. We see little evidence for distributed scattering along most lines of sight in the survey.

 

Terzan 5 is a rich globular cluster within the galactic bulge containing 39 known millisecond pulsars, the largest known population of any globular cluster. These faint pulsars do not have sufficient signal-to-noise ratio (S/N) to measure reliable flux density or polarization information from individual observations in general. We combined over 5.2 days of archival data, at 1500 and 2000 MHz, taken with the Green Bank Telescope over the past 12 years. We created high-S/N profiles for 32 of the pulsars and determined precise rotation measures (RMs) for 28. We used the RMs, pulsar positions, and dispersion measures to map the projected parallel component of the Galactic magnetic field toward the cluster. The 〈B_∣∣〉 shows a rough gradient of ∼6 nG arcsec⁻¹(∼160 nG pc⁻¹) or, fractionally, a change of ∼20% in the R.A. direction across the cluster, implying Galactic magnetic field variability at sub-parsec scales. We also measured average flux densitiesS_νfor the pulsars, ranging from ∼10μJy to ∼2 mJy, and an average spectral indexα= −1.35, whereS_ν∝ν^α. This spectral index is flatter than most known pulsars, likely a selection effect due to the high frequencies used in pulsar searches to mitigate dispersion and scattering. We used flux densities from each observation to constrain the scintillation properties toward the cluster, finding strong refractive modulation on timescales of months. The inferred pulsar luminosity function is roughly power law, with slope$(d\log N)/(d\log L)=-1$at the high-luminosity end. At the low-luminosity end, there are incompleteness effects, implying that Terzan 5 contains many more pulsars.

 

Recently we found compelling evidence for a gravitational-wave background with Hellings and Downs (HD) correlations in our 15 yr data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes, which produce different interpulsar correlations. In this work, we search the NANOGrav 15 yr data set for evidence of a gravitational-wave background with quadrupolar HD and scalar-transverse (ST) correlations. We find that HD correlations are the best fit to the data and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors ∼2 when comparing HD to ST correlations, and ∼1 for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis.

 

Abstract Of the more than 3000 radio pulsars currently known, only ∼300 are in binary systems, and only five of these consist of young pulsars with massive nondegenerate companions. We present the discovery and initial timing, accomplished using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope, of the sixth such binary pulsar, PSR J2108+4516, a 0.577 s radio pulsar in a 269 day orbit of eccentricity 0.09 with a companion of minimum mass 11 M ⊙ . Notably, the pulsar undergoes periods of substantial eclipse, disappearing from the CHIME 400–800 MHz observing band for a large fraction of its orbit, and displays significant dispersion measure and scattering variations throughout its orbit, pointing to the possibility of a circumstellar disk or very dense stellar wind associated with the companion star. Subarcsecond resolution imaging with the Karl G. Jansky Very Large Array unambiguously demonstrates that the companion is a bright, V ≃ 11 OBe star, EM* UHA 138, located at a distance of 3.26(14) kpc. Archival optical observations of EM* UHA 138 approximately suggest a companion mass ranging from 17.5 M ⊙ < M c < 23 M ⊙ , in turn constraining the orbital inclination angle to 50.°3 ≲ i ≲ 58.°3. With further multiwavelength follow-up, PSR J2108+4516 promises to serve as another rare laboratory for the exploration of companion winds, circumstellar disks, and short-term evolution through extended-body orbital dynamics. 

Analyses of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nanohertz frequency band. The most plausible source of this background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for this background and assess its significance make several simplifying assumptions, namely (i) Gaussianity, (ii) isotropy, and most often, (iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated data sets. The data-set length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15 yr data set. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated data sets, even though their fundamental assumptions are not strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them.